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Abstract. - We present compelling evidence supporting the conjecture that the origin of the roton 
in Bose-condensed systems arises from strong correlations between the constituent particles. By 
studying the two dimensional bosonic dipole systems a paradigm, we find that classical molecular 
dynamics (MD) simulations provide a faithful representation of the dispersion relation for a low- 
temperature quantum system. The MD simulations allow one to examine the effect of coupling 
strength on the formation of the roton minimum and to demonstrate that it is always generated 
at a sufficiently high enough coupling. Moreover, the classical images of the roton-roton, roton- 
maxon, etc. states also appear in the MD simulation spectra as a consequence of the strong 
coupling. 
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i\ Since the appearance of the pioneering papers by Lan- 
dau on the theory of superfluidity pQ, the notion of the 
''roton minimum" in the collective mode dispersion of the 
system has played a pivotal role in the explanation of the 
behavior of the superfluid phase of 4 He. This Letter ad- 
dresses this old issue from a novel point of view, made 
possible by results from newly available computer simula- 
tion techniques. 

The revealing feature of the superfluid behavior is the 
non-viscous flow of the liquid (or the frictionless mo- 
tion of an impurity in the liquid [3]) below the critical 
temperature T c = 2.17 K, and the breakdown of this be- 
havior above a critical flow velocity. In Landau's theory, it 
is postulated that in order to explain this breakdown, the 
u(k) dispersion curve of the collective excitation must be 
non-monotonic, starting with a longitudinal u> = ks acous- 
tic portion, reaching a maximum (the "maxon"), which is 
then followed by a deep minimum around a k = fc m i n (the 
"roton minimum"), such that the u = fcs max line, repre- 
senting the relative motion with respect to the fluid, be- 
comes tangent to the dispersion curve in the neighborhood 
of the minimum. 

Since 1947 the physical origin of the roton minimum and 
the question whether it s a part of thelongitudinal phoiion 
eexcitation has been open to conflicting interpretations. 



Landau's original suggestion was that the development of 
the roton minimum had to be sought in the existence of a 
new collective excitation, quite distinct in nature from the 
longitudinal acoustic phonon and rather due to a rotation- 
like collective motion of the fluid: hence, of course, the 
term by which this phenomenon has become to be known. 

Almost a decade later Feynman and Cohen [4] took up 
the detailed analysis of the ground state of a strongly cor- 
related bose liquid. Their seminal work led to the identi- 
fication of S(k), the static structure function, as the cen- 
trally important quantity in the characterization of the 
system. Through a series of physical arguments - whose 
outcome today would be summarized as the "Feynman 
Ansatz" - Feynman was able to relate the ground state 
S(k) to the ground state excitation spectrum ui(k) of the 
fluid through his famous formula 

hk 2 
2mb(k) 

Based on ([1]), the observation that strong correlations cre- 
ate a sharp maximum in S(k) resulted in the re-emergence 
of the roton minimum from a rather different physical 
foundation than that envisioned by Landau. Nevertheless, 
perhaps surprisingly, the interpretation offered by Feyn- 
man still endorsed the association of the roton minimum 
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with a vortex-like excitation. 

In the subsequent decades while observational data 
through neutron scattering measurements accumulated 
and unambiguously revealed the existence of the roton 
minimum at around k m \ n = 1.93 A -1 almost as pre- 
dicted by Landau's hypothesis, the interpretation of its 
physical origin remained somewhat in a limbo. The idea 
that the low-A: phonon and the high-fc roton stem from dif- 
ferent sources prevailed: in a scenario developed by Glyde 
and Griffin in f 990 the roton was viewed as a single parti- 
cle excitation mixed with the density fluctuation spectrum 
[BJ. At the same time, however, in the 1970-s Schneider 
and co-workers [7] pointed at a different role that the ro- 
ton minimum seemed to play: namely being the heralded 
of an incipient liquid-solid phase transition. Still, it took 
an additional thirty years for the corollary of this state- 
ment to be enunciated by Nozieres [5] , by suggesting that 
the physical mechanism responsible for the very creation 
of the roton minimum is the strong correlations prevailing 
in the system (although the fact that local minima ap- 
pear in the dispersion curves of other strongly correlated 
systems resembling the 4 He roton minimum was realized 
much before that [5]). The purpose of this Letter is to 
provide evidence, based on recent computer simulations, 
that favors this latter physical picture. 

New developments in condensed matter physics over 
the past decade, notably the emergence of experimentally 
available novel 3D and 2D bosonic systems, have created 
the feasibility of studying the physics of the roton mini- 
mum under more controlled conditions. Perhaps the most 
interesting new system that has appeared along these lines 
is the two-dimensional liquid of excitonic dipoles. The ex- 
istence of an excitonic liquid in 2D semiconductor layers 
and its expected condensation at T = was predicted 
by Keldysh and co-workers a long time ago [TU|. Subse- 
quently, the possibility of realization of stable excitons in 
electron-hole bilayers was analyzed by Lozovik and others 
[TTj. It is now fairly well established that with a suffi- 
ciently small layer separation d (d/a < 1, where a is the 
average intcrparticlc distance within a layer) the descrip- 
tion of these excitons by a model of a single 2D layer of 
point dipoles interacting through a classical dipolc-dipole 
potential is a very good approximation. 

Recently, the properties of such a 2D system of bosonic 
dipoles, both at zero and at finite temperatures have been 
the subject of intensive analytic and computer simulation 
efforts. Some experimental observations in GaAs semi- 
conductor bilayers |12j can also be interpreted in terms of 
this model. Contrasting with liquid 4 He, the unique fea- 
ture of the dipole system resides in the simplicity of the 
intcrparticle interaction and the possibility of tuning the 
coupling strength by changing system parameters, such 
as temperature and density. It is now feasible to trace 
both the equilibrium behavior and the structure of col- 
lective excitations over a wide range coupling strengths 
and temperatures, by applying various techniques of com- 
puter simulation [T3"1[T4"] . This includes the analysis of the 



generation of the roton minimum and of its evolution as 
a function of the system parameters. These simulation 
studies have shed a new light on this old issue. In the fol- 
lowing we discuss what can be gleaned from these recent 
results and point out what conclusion can be drawn from 
them as to the physical origin of the roton minimum. 

The 2D point-dipole system can be described as a col- 
lection of N spinless point dipoles, each of mass m = m e + 
TTih, occupying the large but bounded area A; n = N/A 
is the average density. The dipoles are free to move in 
the x — y plane with dipolar moment oriented in the z 
direction. The interaction potential is accordingly given 

by 

<p(r)=H 2 /r 3 , (2) 

where fi is the electric dipole strength. The coupling 
strength can be characterized at arbitrary temperature by 

f D = (i 2 /a 3 E kin , (3) 

where £kin is a characteristic kinetic (a = 1/y/im is the 
average inter-particle distance) . At zero temperature, 

E kin ~ h/ma 2 , (4) 

and Td = td = 7*o/a, ?*o = my? /h 2 is the appropriate 
measure of the coupling strength; the characteristic length 
7*o is the dipole equivalent of the Bohr radius. In fact, (0| 
is valid up to an arbitrary factor 0(1): its exact evalu- 
ation would require information about the ground state 
momentum distribution. In the high-temperature classi- 
cal domain, with 

-Ekin = (5) 

(/3 _1 = k-oT is the thermal energy) and ^ becomes 
Td = fll^/a, 3 , in analogy with the T = f3e 2 /a coupling 
parameter for classical Coulomb liquids. The comparison 
of the two limits of Td allows one to establish the cor- 
respondence between the high temperature classical) and 
low temperature (quantum) domains through the equiva- 
lence Td I'd- 

The behavior of the dipole system can be analyzed over 
a wide range of temperatures and densities. The temper- 
ature domains may be identified through the value of the 
degeneracy parameter 

6 = f3h 2 /2ma 2 . (6) 

As the quantum, > I and classical, < I and 
regimes. The equilibrium in the quantum domain ex- 
hibits a great richness of phenomena, such as condensation 
and quasi-condensation, superfluidity and crystallization, 
whose competition leads to an intricate low temperature 
phase diagram [15], quite at variance with the behavior 
in the classical domain. However, as far as the collective 
excitation spectrum is concerned, the situation is differ- 
ent: collective modes of many-body systems governed by 
the cooperative motion of many particles seem to be ro- 
bustly insensitive to statistics (an obvious example is the 
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plasmon in charged particle systems, which exhibits very 
much the same features in classical plasmas as in the de- 
generate electron gas). 

We proceed now to compare the description of the col- 
lective mode dispersion based on two sets of computer 
simulation data. To the first set belong data pertaining to 
the static equilibrium properties of low temperature quan- 
tum dipole systems, through zero temperature quantum 
Monte Carlo (QMC) [T3] and low-temperature path inte- 
gral Monte Carlo (PIMC) [T5] simulations. The second 
set of data originates from a Molecular Dynamics (MD) 
simulation of the dynamics of a classical point dipole sys- 
tem. The QMC simulation of the zero temperature quan- 
tum dipole system in the strong coupling regime has been 
pursued by Astrakharchik ct al. (T3] and by Mazzanti et 
al. [16]. In these works static structure factors have been 
generated as functions of the zero temperature coupling 
parameter r^. The dispersion relation for the longitudi- 
nal dipole oscillations was obtained through two analytic 
approximations: first, through the Feynman formula ((1}, 
and, second, by incorporating the Feynman frequencies 
in a more sophisticated formalism based on the correlated 
basis function (CBF) approximation (extensively reviewed 
by Feenberg [T7], Campbell [TS], and Woo [H5]), which 
takes account of the effect of the three-phonon interaction 
term in the Hamiltonian. In the PIMC simulation ap- 
proach [20] followed by Filinov and co-workers [15] , similar 
data for low but finite temperatures have been generated. 
Coming now to the second set of data, a program of MD 
simulation of the dynamics of a classical point dipole sys- 
tem has been pursued by the authors, along the pattern 
used for simulating other strongly interacting classical sys- 
tems |21l . In this approach the dynamical structure func- 
tion S(k,uj) is obtained for a range of values of the clas- 
sical coupling parameter [22j To whose peaks reveal the 
collective excitation spectrum of the system. It should be 
emphasized that, in contrast to the work of Astrakarchik 
et. al. |13j . this latter method provides direct informa- 
tion without the intervention of any further analytic ap- 
proximation on the dispersion relation for the longitudinal 
dipole oscillations in the physical, albeit classical, system. 
Within the presently available computational technology 
it would be of course futile to expect data of the same kind 
to be created for quantum systems, where, to date, one has 
to rely on the intermediary of approximation techniques 
to connect the dynamics of the system with the available 
equilibrium information. 

We now compare the results from the two sets of data, 
which are presented in Fig. [TJ for approximately equiva- 
lent values of the classical and zero temperature coupling 
parameters T d . The main statement of this paper derives 
from a number of conclusions that can be drawn from the 
inspection and analysis of these graphs. 

1. There is no major difference in the morphologies of 
the dispersion curves pertaining to the quantum and 
classical systems; in particular, both of them show 
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Fig. 1: Comparison of dispersion curves: MD (classical molecu- 
lar dynamics) , Astrakarchik [13] (Feynman upper bound based 
on QMC), Mazzanti [16] (presumed lower bound), (a) Yd = 
15, (b) T D = 30. 

the formation of the roton minimum in the vicinity 
of ka ~ 3.5. 

2. Whether the roton minimum develops in the disper- 
sion curve depends solely on the coupling parameter 
r o : the minimum value for the curve to become non- 
monotonic is around = 5, both for the classical 
and quantum system. 

3. With increasing value of T^, the roton minimum be- 
comes deeper, reaching an absolute minimum. For 
the T = curve this occurs at Yd = 28.6, near the 
crystallization value [13] Y* D = 30.2. 

4. The classical dispersion curve is sandwiched between 
the quantum Feynman and the quantum CBF curves; 
its being positioned below the Feynman curve is ex- 
pected, since the latter constitutes an upper bound 
only. As to the CBF construction, it is difficult to 
assess its accuracy by analytic means; however, CBF 
calculations performed for 4 He with model potentials 
have agreed well with measured spectra, with devi- 
ations attributable to the uncertainty of the actual 
interaction potential [23j . Thus, one is probably on 
safe grounds by accepting the CBF curve as a lower 
bound. 

All these observations converge now towards a coherent 
physical picture that clearly suggests that the roton min- 
imum is the consequence of strong correlations, which are 
basically a classical phenomenon, in the system. Quantum 
dynamics, Bose statistics, condensation, superfluidity, etc. 
seem to have very little influence either on the formation 
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Fig. 2: Comparison of T = and finite low temperature dis- 
persions. Note that the finite temperature modifications are 
very small. 



Fig. 3: Broadening of the phonon-maxon-roton line due to tem- 
perature effect. Error bars indicate the width (full width at half 
maximum) of the spectral peak. 



or on the structure of the roton minimum. It would be a 
strange coincidence indeed, if hybridization with vortex- 
like or other unidentified excitations were responsible for 
the creation of the same phonon dispersion characteris- 
tics in bose liquids as they are manifested through strong 
correlations in the equivalent classical systems. 

Having made our central observation, there arc a few 
more issues to be addressed. Our discussion on the quan- 
tum system so far has been limited to the T = ground 
state, whose analysis is sufficient to establish the correla- 
tional origin of the roton minimum. Finite temperatures 
are, however, now available: the recent PIMC calculations 
[15] have generated S(k) data in the 0.5 < 9 < 3.3 do- 
main (on both sides of the superfiuid/normal fluid tran- 
sition point). The generalization of the zero temperature 
Feynman formula ([T} to finite temperatures becomes 



)(k) tanh 



/?7kj(fc) 



fik 2 



2mS(k) 



(7) 



The PIMC dispersion curves based on ([7]) for three rep- 
resentative values, together with the QMC ground state 
dispersion are displayed in Fig. [5J The difference between 
the T = and T ^ dispersions is strikingly small (ex- 
cept near k — 0, a region for which the inaccuracy of the 
QMC simulations [13] has already been pointed out [21] . 
further corroborating the expectation that as far as the 
dispersion properties of the collective excitations are con- 
cerned, the transition from the T = degenerate to the 
T — > oo classical phase is smooth and eventless. 

Quite at variance with the dispersion, one expects that 
the lifetime of the excitations (i.e. the damping of the os- 
cillations) is strongly affected by rising temperature. Here 
one has to rely on what is known from observations on liq- 
uid 4 He, which then can be compared with relevant infor- 
mation gleaned from the classical MD simulations. From 
the analysis of S(k, uj) with the aid of neutron scattering 
experiments [25] we know that the character of the roton 
minimum peak undergoes a dramatic change as the system 
passes from the superfluid to the normal phase: while it is 
very sharp in the superfluid phase, it becomes broad in the 
latter and as the temperature continues to rise above its 



critical value it becomes quite ill-defined. Thus one would 
expect a similarly broadened roton peak in the classical 
MD generated S(k,uj). This indeed is the case. An ex- 
amination of Fig. [3] that portrays S(k,tu) and the width 
of the roton peak for a sequence of increasing Tp values 
shows that the width decreases for higher coupling, yet 
remains considerable up to crystallization. Beyond ob- 
serving the qualitatively identical behavior in the classical 
2D dipole system and 4 Hc, quantitative comparison would 
be quite meaningless, first because of the entirely different 
interaction potentials in the two systems, and, second, be- 
cause studying the effect of the variation of the coupling 
strength is obviously not an available option in 4He exper- 
iments. The conclusion to be drawn from the comparison 
is, however, that the temperature engendered damping of 
the rotons is similar in the classical and quantum systems 
and its occurrence in no way invalidates the link between 
the finite T classical and the T ~ quantum systems. 

If we now accept the premise that the roton minimum 
experienced in the bose systems studied are caused by 
the strong correlations of classical origin, one may won- 
der about the generality of the link between strong cor- 
relations and roton-like behavior for other many- body 
systems. First, would systems governed by interaction po- 
tentials other than the ones that have been investigated 
exhibit a similar behavior? With a high level of confi- 
dence one can predict that this would indeed be the case. 
In fact, it was prior to the MD analysis of the dipole sys- 
tem that the classical " roton minimum" phenomenon was 
established for Coulomb and Yukawa systems [3T]. There 
is no reason to believe that the mapping of this behavior to 
the equivalent low temperature Bose system, should such 
a system become available for observation, would not oc- 
cur in the same fashion as it does for a system of dipoles. 
Second, could the classical roton minimum behavior be 
mapped to a strongly interacting low temperature system 
of fcrmions in the same way as it is to a system of bosons? 
We do not have enough observational or theoretical un- 
derstanding of such systems to provide a definite answer. 
Nevertheless, it is clear that the fundamentally different 
low temperature behaviors of the fermi and bose systems 
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would make a similarity unlikely. The issue is that in 
a bose condensate or quasi-condensate the single particle 
excitations are suppressed, which makes the unimpeded 
development of the collective excitation, even in the high- 
k domain where the roton minimum is to occur, possible. 
There is no such suppression in the fermi system. Tak- 
ing the electron liquid as a paradigm, the high-fc devel- 
opment of the plasmon branch is profoundly affected by 
the presence of the electron-hole pair continuum. Based 
on a simplified model |26j , with a measure of caveat it can 
be argued that the pair (and multi-pair) continuum not 
only generates a strong damping to the collective mode, 
but also inhibits its evolution in the high-fc domain. Thus 
it seems safe to conclude that the condensed (or quasi- 
condensed) bose system is the sole candidate for the clear 
manifestation of the correlation induced roton minimum. 

The final issue we want to address is what has become 
to be known in the 4 Hc studies as the issue of the two- 
roton state. It was shown first by Cowley and Woods [2"T] 
that in addition to the phonon excitation, a high frequency 
branch can be identified in the observed mode spectrum 
of the liquid 4 He: this mode was identified as a two roton 
excitation, with a frequency in the vicinity of twice the 
frequency of the roton minimum and rather insensitive to 
the wave number. More recent neutron scattering studies 
|28j have found additional branches with higher frequen- 
cies that can be attributed to roton-maxon and maxon- 
maxon excitations. 

The identification of these high frequency states in 
terms of the excitation of two quasi-particlcs goes back to 
the classic paper of Miller, Pines and Nozieres [29]. In the 
original model, the quasi-particles were non-interacting: 
the effect of the interaction between them was suggested 
by Pitaevski [3U] and first considered in detail in the pi- 
oneering series of works by Ruvalds, Zawadowski, Bedell 
, Pines and collaborators [30] • In these works, depending 
on the sign of the interaction (attractive vs. repulsive), a 
downward or upward shift from the simple additive nom- 
inal values of the high frequency excitations is predicted. 

The generation of higher harmonics and combination 
frequencies is a well-known nonlinear effect [31] . due to 
the inherent nonlinearity of the interaction between the 
constituent ( 4 He, dipoles e.g.) particles. Given a disper- 
sion curve, any pair of frequencies along the curve is a 
candidate for the generation of combination frequencies: 
it is, however, the frequencies in the vicinity of extrema 
of the dispersion curve, where the density of states is the 
highest, which are primarily responsible for creating ob- 
servable new branches of collective excitations. In this pa- 
per we have argued that the generation of the roton mini- 
mum and of the maxon maximum in strongly coupled bose 
condensate (quasi-condensate) liquids is a universal ef- 
fect due the strong interaction between their constituents. 
A corollary to this statement is now that the roton- ro- 
ton, roton-maxon and maxon-maxon states should also be 
identifiable for the equivalent classical liquid through the 
appearance of a significant weight in the density correla- 




Fig. 4: (color online) Illustration of the appearance of the 
R+R, R+M, M+M combination frequencies in the longitudinal 
current-current correlation functions L(k,uj), measured along 
the direction of the nearest neighbor in a hexagonal lattice). 
The center color map of the L(k, ui) fluctuation spectra shows 
the strong primary dispersion and the "ghosts" of the combina- 
tion frequencies. The small panels show vertical cross-sections 
taken for the selected group of wave-numbers, where the ap- 
pearance of the combination frequencies is the most manifest. 



tion spectra at the corresponding combination frequencies. 

The MD analysis of S(k,ui) of the strongly coupled 
dipole liquid when extended into the high frequency 
regime up to and beyond twice the maxon frequency in- 
deed reveals these expected features. In Fig. 2] we have 
displayed a series of S(k,ui) spectra for a wide range of 
k values. In order to emulate the low level of the back- 
ground noise characteristic of the superfluid, the graphs 
were taken at To = 500, well beyond the crystallization 
value, but at a value of coupling where the thermal mo- 
tion is of sufficient amplitude to generate nonlinear effects. 
Thus these graphs may be considered appropriate to illus- 
trate a trend, but a detailed analysis of the features exhib- 
ited in them is not our purpose here: it is sufficient to point 
out that the essential feature, namely the accumulation of 
weight in the vicinity of three combination frequencies, 
the roton-roton (R+R), roton-maxon (R+M), and maxon- 
maxon (M+M) is clearly visible. From other works we 
know that the generation of harmonics, similarly to the 
development of the roton minimum itself, appears to be a 
general feature of strongly coupled liquids with Yukawa, 
Coulomb and other similar types of interaction |32j as well. 
The relative amplitudes of these harmonics are very much 
functions of the particular type of interaction potential 
prevailing in the system. Thus, attempts at a quantitative 
comparison with experimental data on liquid 4 He would in 
all likelihood not be useful. On the other hand, the posi- 
tions of the combination frequencies, i. e. their shifts from 
their nominal values, is due, in the language of quasipar- 
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ticle states, to the interaction between the quasiparticles 
and, as such, their comparison with the predictions of the 
RZ theory as applied to the 2DDS would certainly be of 
interest. 

To summarize, in this paper we have provided com- 
pelling support for two earlier suggestions [7J[5], that 
the formation of maxon-roton excitations in condensed 
or quasi-condensed bose systems is the result of a basi- 
cally classical strong correlational effect. Our assessment 
is based on the analysis of the collective mode dispersion 
of the strongly coupled 2D bosonic dipole system (2DDS) 
below and above the superfluid transition temperature 
13, 14]. This dispersion exhibits a phonon-maxon-roton 
structure, very much the same as the collective mode dis- 
persion in 4 He. Paralleling these studies, we investigated 
[2"2"] the dispersion of the strongly coupled 2DDS in the 
high-temperature classical domain through molecular dy- 
namics (MD) simulations: we have discovered a remark- 
able congruence between the collective mode dispersions of 
the classical and quantum systems, pointing at a common 
origin of these phenomena, the shaping of the longitudi- 
nal phonon dispersion by strong intcrparticlc correlations. 
Models that seek the origin of the roton in excitations dif- 
ferent in nature from the phonon seem to be incompatible 
with this evidence. Moreover, the MD simulations now 
can be trusted to provide a first-time glimpse into the ac- 
tual behavior of the dynamics of the strongly correlated 
system without reliance on the intermediary of models. 
The approach that describes correlational effects in the 
dynamical behavior through the static structure function 
is superseded by the direct dynamical MD simulations. 
Furthermore, our observations of a considerable weight in 
the MD density fluctuation spectrum in the vicinity of the 
roton-roton, roton-maxon, and maxon-maxon combina- 
tion frequencies, similar to what has been observed in neu- 
tron scattering experiments [55] in 4 He, may be regarded 
as additional corroboration of the underlying strong cor- 
relational classical model. 
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